Magneto-oscillations due to electron-electron interactions in the ac conductivity of a 

2D electron gas 
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Electron-electron interactions give rise to the correction, Sa lnt (uj), to the ac magnetoconductivity, 
o~(lu), of a clean 2D electron gas that is periodic in to' 1 , where lo c is the cyclotron frequency. Unlike 
conventional harmonics of the cyclotron resonance, which are periodic with u>, this correction is 
periodic with oj 3//2 . Oscillations in 8a lnt (ui) develop at low magnetic fields, ui c <C u, when the 
conventional harmonics are suppressed by the disorder. Their origin is a double backscattering of 
an electron from the impurity-induced Friedel oscillations. During the time ~ between the two 
backscattering events the electron travels only a small portion of the Larmour circle. 

PACS numbers: 73.40.-c, 73.43.-f, 73.43.Qt, 78.67.-n 
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Introduction. Originally, the cyclotron resonance (and 
its harmonics) in the ac conductivity, <j(uj), of the 2D 
electron gas had been detected by measuring the trans- 
mission of the microwave radiation wi In the recent 
experiment on high- mobility samples |2], it was demon- 
strated that this resonance, together with harmonics, also 
manifests itself in the dc magnetoresistance under mi- 
crowave illumination, i.e., in the photoconductivity. A 
spectacular strength of this effect, and, in particular, ob- 
servation of the zero-resistance state, above a certain in- 
tensity of illumination 0, Q, had attracted a steady 
interest of the researchers to the ac-response of a high- 
mobility electron gas in a weak magnetic field, B. Unlike 
the conventional Shubnikov-de Haas oscillations of the 
dc magnetoresistance, which vanish with temperature as 
cxp(— 2n 2 T / uo c ), where uj c is the cyclotron quantum, the 
magneto-oscillations of <j(uj) survive at high temperature 
d, 0] • The shape of these oscillations is given by 0, 0| 



5<t±(u) „ / 2trj\ / 2tt 
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where r is the scattering time, and a± (ui) is related to the 
dc conductivity <tq as a± (w) = a /[2 + 2{u±iu c ) 2 T 2 ]. 
Classically, the meaning of the damping factor S 2 — 
exp(— 2n/uj c T^ < 1 is the probability for an electron to 
execute the entire Larmour circle, 27ri? L , without being 
scattered. Oscillations Eq. (p} is a single- electron effect. 
In converting of these oscillations into the oscillating dc 
photoconductivity @, H, H, OH , the electron-electron in- 
teractions enter as a source of relaxation of the oscillatory 
part of the distribution function. 

In the present paper we demonstrate that interactions 
by themselves give rise to the oscillatory contribution, 
5a lnt (u), to the linear ac conductivity, cr(u)), at frequen- 
cies much higher than in Eq. ([I]). To contrast this con- 
tribution to Eq. {l}, we present 8a mt {Lo) in the form 
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FIG. 1: Illustration of the S-dependent contribution, SQ B (r), 
to the phase of the polarization operator; 7? L [# — 2sin(6>/2)] 
is the elongation of the semiclassical trajectory due to the 
field-induced curving. The origin of oscillating magnetocon- 
ductivity is the scattering from the Friedel oscillations (arcs 
of decreasing thickness), at points 1 and 2, located symmet- 
rically with respect to the impurity shown with a big dot. 



6a int (tu) ( cj vr\ 
-— oc cos \C U exp 



^( l - + 2,T 



where 



a„ = 
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and E F is the Fermi energy. Since C u is small, the correc- 
tion Eq. ^ develops oscillations at much smaller mag- 
netic fields lo c ~ C^juo <C ui than Eq. ([1]). At such fields, 
the damping factor in Eq. (JTJ) is ~ exp (— 2-k/CJ), i.e., 
the conventional oscillations are completely washed out. 
Qualitative picture. The origin of the oscillations Eq. ^ 
lies in a peculiar modification by the interactions of the 
impurity scatte ring in a weak magnetic field. Conven- 
tionally, ll|, \12l 1 1 311 this modification amounts to the ad- 
ditional scattering [14[ from the Friedel oscillations of the 
electron density, created by the impurity. Such a modifi- 
cation does not lead to the anomalous sensitivity to low 
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FIG. 2: Diagrams contributing to the second-order interaction 
correction to the ac magnetoconductivity. Dashed lines de- 
note the impurity scattering. Dots in the vertices combine two 
types of the interaction matrix elements, as shown in the left 
inset. Diagrams b and e describe the impurity scattering of 
the secondary electron (hole) ; other diagrams describe double 
scattering of photoexcited electron (hole) by the Friedel oscil- 
lation. Right inset: diagram (a) for purely disorder-induced 
part of the self-energy, E n , is plotted together with represen- 
tative diagrams (b) and (c) for the interaction-induced self- 
energy, JEJj' 1 '. Only diagrams (c), describing two electron- 
electron scattering processes, contribute to the interaction- 
induced magneto-oscillations. 

B. However, as we demonstrate below, this sensitivity 
emerges in the second order in the electron-electron in- 
teraction strength. The corresponding second-order pro- 
cesses are illustrated in Figs. 1 and 2. They are: (i). 
Photoexcited electron emits a virtual pair, which is sub- 
sequently annihilated. Impurity scatters not the original 
electron, but rather the impurity scattering occurs be- 
tween the states, constituting the pair, prior to annihila- 
tion. Diagram b in Fig. 2 describes this process, (it) Elec- 
tron is not scattered directly by the impurity, but rather 
experiences a double backscattering from the impurity- 
induced Friedel oscillations, as illustrated in Fig. 1, and 
also by the diagram a in Fig. 2. 

We demonstrate that the corrections to the conduc- 
tivity from both these processes oscillate with magnetic 
field according to Eq. ([2|) . The oscillations reflect the fact 
that, for both processes, the dominant contribution to the 
double backscattering cross-section comes from two "dis- 
tinguished" points that are located symmetrically with 
respect to the impurity at certain well-defined distance, 
r w ~ C^i? L , see Fig. 1. Then the argument of the cosine 
in Eq. ^ can be interpreted as a product wt u , where 
tui = r u/v F is the time, during which the electron with 
Fermi velocity, v F , travels the distance r u . 
Derivation. Although the interaction-induced oscilla- 
tions come from small distances, r w <C i? L , we neverthe- 
less will evaluate o~(uj) in the Landau gauge to demon- 
strate how both oscillations Eq. ([1]) and Eq. ([2]) emerge 
from the same calculation. Within the self-consistent 
Born approximation, averaging in the general expression 
for the diagonal conductivity 



a ( UJ ) = a in, / — (/« ~ /e+w) Tru x ImG e+w u x ImG 6 (4) 
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is decoupled into two averaged Green functions 

ft = £ i (5) 

n n v ' 

where e„ = (n + ^)hjj c are the Landau levels, and £ n (e) 
is the self-energy. In Eq. (01, the bar denotes the dis- 
order averaging, O is the normalization area, and f e is 
the Fermi distribution. Upon decoupling, Eq. takes 
a familiar form 

e 2 vluj c v fde f \ 
x^ImG n±1 (e + W )ImG„(e), (6) 

n 

where i^o is the 2D density of states. For high Lan- 
dau levels, (E F /huj c ) 3> 1, in the first approximation 
in S < 1, the self-energy can be replaced by its zero- 
field value, z/2t. Then Eq. © readily reproduces the 
Drude conductivity. In order to capture the oscilla- 
tory ac magnetoconductivity, in the next approximation, 
one should take into account the "quantum" correction, 
5T,Q(e) oc (5exp (— 2TTie/uj c ), to the self-energy due to the 
discreteness of the Landau levels, as well as the inter- 
action correction, <5£™*(e). Since both corrections are 
smaller than 1/r, they cause a small correction to the 
Green functions Eq. ([5]) of the form 

. «3M+"r . (7) 

(e - e„ - £) 

The first and the second terms in Eq. ([7|) give rise to the 
oscillations Eq. ([I]) and Eq. ([2|), respectively. However, 
to reproduce these oscillations the "quantum" and the 
interaction corrections should be handled differently. To 
reproduce Eq. ([1]), upon substituting Eq. ([7]) into Eq. ([6]), 
one should keep the product, <5£^ +1 (e + oj) [5E$(e)l*. It 
contains the oscillating term oc exp (— 27riw/w c ), which 
does not depend neither on n nor on e. For this rea- 
son, the resulting oscillations of magnetoconductivity are 
T- independent. By contrast, to capture the interaction- 
induced oscillations, it is sufficient to keep <5E" lt only 
in one of the Green functions in Eq. ([S]), and its n- 
dependence is crucial. We will perform further calcu- 
lation for <5E™* (e) given by the first diagram of type c in 
Fig. 2 (inset). This is because the diagrams of type b do 
not cause magneto-oscillations, while the contributions 
of other diagrams of type c are comparable to that of the 
first one, and will be addressed later. 
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The first diagram of type c can be presented as 

6K nt (e)=Yl lRnm [ l , (8) 



so that the n-dependence is encoded in the "matrix ele- 
ments", R nm - Substituting Eq. ([8]) into Eq. (JT]), and then 
Eq. into Eq. © yields 



5a mt {cu) oc 



(9) 



I Rnm I 



n,m (e - e« + ^p) (e + a; - e n ±\ - j^) (e - e m + 27) 

As a next step, we express the matrix element, R nm , as 
an integral in the coordinate space, following Fig. 2c, 
R nm oc / <irV'*(r)'!/'m(r)no(r, 0), where n (r, 0) is the 
static polarization operator between the point r = 0, 
where the impurity is located, and the point r, where 
the backscattering takes place. Our prime observation is 
that with such R nm the relevant term in Eq. ([9]), which 
has the form 



de (f e+u - f e ) 

(e + w - e n±1 -ip)(e-e m + 



(10) 



again reduces to the polarization operator, n w zp We (ri, r 2 ). 

In the final expression for the interaction correction 
we make use of the fact that the -B-dependence of this 
correction develops in the low-field limit oj c ~ C^oj <C oj, 
and replace u> ± cj c by u> 15]. We then obtain 



5a 



X 2 



[dr 1 [dr 2 Imn w (ri,r 2 ) 
xRe {n (0,ri)n (r 2 ,0)}, (11) 

where A stands for dimensionless strength of interaction, 
which we assumed to be short-ranged. The numerical 
factor in Eq. (jTTJ) will be established when all the dia- 
grams contributing to 5a mt are considered (see below). 
Interpretation. The form of Eq. (|lip can be interpreted 
as follows. The factor, ImILj(ri, r 2 ), in the integrand 
is the density-density response, the same as in calcula- 
tion of the Drude ac conductivity. The second factor, 
Re |lTo(0, ri)Ilo(r2, 0)|, plays the role of the spatial cor- 
relator of the effective random potential. By lifting the 
momentum conservation, this potential enables the ab- 
sorption of the ac field. If the correlator was oc 8{yi — r 2 ), 
then the rhs of Eq. (fTTj) would yield an w-indepcndent 
constant. Important is that the effective potential in 
Eq. originates from the modulation of the electron 
density by the impurity, and thus oscillates rapidly with 
distance. It is these Friedel oscillations that in magnetic 
field lead to the oscillating correction, Eq. @. 



Oscillations. The long-distance, k F r 3> 1, behavior of the 
polarization operator in coordinate space is the following 



n w (r) 



' 2k F r 



sm{6(r)} / 



x exp< 



— , (12) 

V F T 



where the function A(x) — x/ smh(x) describes the tem- 
perature damping. In the momentum space, two con- 
tributions to Eq. (|12[) originate from small momentum 
transfer and momentum transfer close to 2k F , respec- 
tively (l6| . At distances r <C i? L , a nonquantizing mag- 
netic field enters into Eq. (TT2"|) through the semiclassi- 
cal phase, 0(r). This phase is accumulated by the elec- 
tron upon propagation from the point to the point r 
and back. In a zero magnetic field, we obviously have, 
0(r) = 2k F r. At distances r <C i? L , the field-dependent 
correction [l7| to 6(r) is equal to 



AB 

6@ B (r)=2k F 5£ =-- 

$ 



6i£ 



(13) 



The origin of the correction Eq. ([13)) is illustrated in 
Fig. 1. It comes from elongation, 5C = R L {9 — 2sin(#/2)], 
of the classical electron trajectory in magnetic field, as 
well as from the Aharonov-Bohm flux into the loop with 
area„4 = (9 — s'm9)R^/2. The correction Eq. (fT3"|) is nega- 
tive, since the Aharonov-Bohm contribution exceeds twice 
the orbital contribution. We emphasize, that the conven- 
tional way [18] of incorporating magnetic field into the 
Green's function neglects the curvature of the electron 
trajectories, i.e., 5<d B (r) = 0. Thus, within the approach 
of Ref. [lg, the oscillations Eq. @ would not emerge. 

Further calculation is straightforward. Substituting 
Eq. (fT2|) into Eq. (jTTJ) , performing the angular integra- 
tion, and combining rapidly oscillating terms in the prod- 
uct of three polarization operators into a "slow" term, 
we find that the interaction correction Eq. (TlTj) can be 
presented as 5o- mt /a(u}) = (A 2 E F /u)F TjT , where the di- 
mensionless function F Tj t(cj, lo c ) is defined as follows 



dr\dr 2 



(7rfc P ) 5/2 J hr 2 (r 1+ r 2 )] 3/2 V « 



A {^) A { Mri v F r2)T )^{ 
LL>(ri+r 2 ) 



2{n+r 2 ) 



x(14) 



■ + S6 B (r ! + r 2 ) - <5© B (r 1 ) - <5© B (r 2 ) - 



Other slow terms emerging in the rhs of Eq. (jTTJ) , 
e.g., the one with uj — » — lu, do not oscillate with mag- 
netic field. By contrast, the function F T .r does oscil- 
late, since the argument of cosine in Eq. (I14p . with 5<d B 
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FIG. 3: (Color online) Interaction-induced contribution to the 
ac conductivity calculated numerically from Eq. (|14p (dotted 
line) and from asymptotic expression Eq. (|16p (full line) are 
plotted vs dimensionless frequency x = 2 1 ^ 3 'lo / '(E F ^ 3 'uiV^) = 
3(z/4) 2 / 3 . The calculations are performed for dimensionless 
disorder 1/(E f t) = 0.08(tJ c /i?F) 2/ ' 3 and dimensionless tem- 
perature T/E F = 0.06(lo c /E f ) 2/3 . 



given by Eq. (|13p. has a saddle point at r± = r 2 = r u = 
3C u v F /4w c = (3/4) C U R L . In the vicinity of the saddle 
point, the phase of the cosine can be presented as 
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(ri -r,) 2 + (r 2 - n,) 2 
. (ri - r u )(r 2 - r u ) 



(15) 



where is defined by Eq. (jSJ). Note, that the combi- 
nation, [C^cj/oJc — f] , in (|15p is nothing but the phase 
of the interaction-induced oscillations Eq. @. It also 
follows from Eq. (fl~5|) that, when this phase is large, the 
characteristic deviations, (r% — r w ) and (r 2 — r w ) are much 
smaller than r w . This allows to perform the integration 
over these deviations in Eq. (I14p explicitly. This yields 



F TjT (cj,cj c ) = 




(3 • 2 29 / 3 7T 3 ) 



cos( 9^L - 1 )a 2 ( 2?rr " r \l f 4?rr " r 
\ w c 4 / \ v F J \ v F 



exp 



4r^ 



Other diagrams. Contribution Eq. (f!6| to (5er m *(tj) is 
the result of calculation of a single diagram a in Fig. 2. 
Other diagrams, involving two electron-electron scatter- 
ing processes and yielding contributions with a structure 
similar to Eq. (fl7J|) , are shown in Fig. 2. Diagrams b, c, 
and d are captured within the self-consistent Born ap- 
proximation, and correspond to certain terms in ST, mt , 
see c) in Fig. 2 (inset). Diagrams e-h in Fig. 2 are of 
the same order as a-d, but they are not contained in 
Eq. (|6j); these diagrams emerge from the general expres- 
sion Eq. Q for cr(u>). Taking all the diagrams into ac- 
count leads to the modification of the function D(z) from 
z~ 5 / 6 to D(z) = -32.Z- 5 / 6 + 64Z 1 / 6 , where the factors 
—32 and 64 account for the spin indices and for the num- 
ber of closed fermion loops in different diagrams. The 
second term, oc z 1 / 6 , arises from the diagrams b, e and 
d, h in Fig. 2. Since the oscillations in Eq. ([2]) develop at 
z»l, these diagrams are, actually, dominant. 

Numerical estimates. Note that, in terms of B- 
periodicity, oscillations Eq. ([2]) coincide with oscillations 
Eq. ([1]) upon rescaling oj c by 2tt/C m in the argument of 
cosine, and by 27r/3C w in the Dingle factor. For a typi- 
cal ac frequency hw ~ 3K and density n ~ 10 n cm _2 in 
the experiments 0, 0, 0, 0, EH, IS] this shifts the domain 
of oscillations Eq. © from B ~ 0.2T to B < 1(T 2 T. 
For such B the observation of the oscillations requires 
T < lUc/QttCu ~ 20mK, which was not the case in 
Refs. [3, i, i, H, [H, 0. For observation of magneto- 
oscillations Eq. ([2]) higher densities n ~ 5 ■ 10 n cm~ 2 and 
frequencies fiw ~ 15K are needed. 



where D(z) = z~ 5 / 6 . For high enough temperatures 
T > 2uJc/3ttCuj (but still T<u), the damping factor can 
be replaced by the exponent, and we reproduce the os- 
cillating contribution Eq. ([2]). The above derivation sug- 
gests that oscillatory behavior of the correction Eq. @ 
establishes only at large z 3> 1, which corresponds to 
the nodes of cos(z — 7r /4) with high numbers. To find 
out where the asymptotics Eq. (fT6|) actually applies, we 
have evaluated the double integral Eq. (fT4|) numerically. 
The result is plotted in Fig. 3 and indicates that Eq. @ 
applies starting already from the third node. 
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